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We review our recent study of the pion mass dependence of the p and a resonances, generated 
from one-loop SU(2) Chiral Perturbation Theory (ChPT) with the Inverse Amplitude Method 
(IAM) which was modified to properly account for the Adler zero. The method is based on 
analyticity, elastic unitarity and ChPT at low energies, thus yielding the pion mass dependence 
of the resonance pole positions from the ChPT series up to a given order. We find that the pitit 
coupling constant is almost m K independent and that our prediction compare well with some 
recent lattice results for the p mass. These findings may be relevant for studies of the meson 
spectrum and form factors on the lattice. 
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1. Introduction 

Light hadron spectroscopy at low energies lies beyond the realm of perturbative QCD. Al- 
though lattice QCD provides, in principle, a rigorous way to extract non-perturbative quantities 
from QCD, current lattice results are typically still obtained for relatively high quark masses. Thus, 
in order to make contact with experiment, appropriate extrapolation formulas need to be derived. 
This is typically done by using Chiral Perturbation Theory (ChPT)[j|], which provides a model 
independent description of the dynamics of the lightest mesons, namely, the pions, kaons and etas, 
which are identified with the Goldstone Bosons (GB) associated to the QCD spontaneous Chiral 
Symmetry breaking. Hence, ChPT is built out of only those fields, as a low energy expansion of a 
Lagrangian whose terms respect all QCD symmetries, and in particular its symmetry breaking pat- 
tern. Actually, this chiral expansion becomes a series in momenta and meson masses, generically 
0(p 2 /A 2 ), when taking into account systematically the small quark masses of the three lightest 
flavors that can be treated perturbatively. The chiral expansion scale is A = Anf n , where f % de- 
notes the pion decay constant. ChPT is renormalized order by order by absorbing loop divergences 
in the renormalization of parameters of higher order counterterms, known as low energy constants 
(LEC), which are the coefficients of the energy and mass expansion, so that they have no quark 
mass dependence. Their values depend on the specific QCD dynamics, and have to be determined 
either from experiment or from lattice QCD. 

The relevant remark for us is that, thanks to the fact that ChPT has the same symmetries as 
QCD and that it couples to different kind of currents through the same pattern, the ChPT expansion 
provides a systematic and model independent description of how the observables depend on some 
QCD parameters, like the quark masses, and this can be implemented systematically up to the 
desired order in the ChPT expansion. 

We review here our recent derivation of a modified version of the IAM [Q] based on dispersion 
theory, unitarity and ChPT to next-to-leading order (NLO), which also accounts properly for the 
Adler zero. Within this approach we are able to predict the quark mass dependence of the a and p 
mesons [g], thus providing an explicit representation of the LECs appearing in the ChPT analysis 
of the vector meson mass of Refs. [f|, g]. 

In this work we focus on the two lightest resonances of QCD, the p and the a. It is therefore 
enough to work with the two lightest quark flavors u,d in the isospin limit with a mass m = (m u + 
md)/2. Since m n is given by m\ ~ m+ ... [|J], studying the m dependence is equivalent to study the 
m % dependence. 

2. Unitarization and dispersion theory 

The a and p resonances appear as poles on the second Riemann sheet of the (/,/) = (0,0) 
and (1,1) Tin scattering partial waves of isospin / and angular momentum /, respectively. Unitarity 
implies for these partial waves, and physical values of s below inelastic thresholds, that 

lmt(s) = a(s)\t(s)\ 2 =>- Im— = -ct(j), with a(s)=2p/y/s, (2.1) 

t(s) 

where s is the Mandelstam variable and p is the center of mass momentum. Consequently, the 
imaginary part of the inverse amplitude is known exactly. Hoewever, ChPT ampitudes, being an 



2 



Chiral Extrapolations of light resonances from dispersion relations and Chiral Perturbation Theory 
Guillermo RIOS MARQUEZ 



expansion t ~ t 2 + k H , with tk = 0(p 1 ), can only satisfy Eq. (2. 1 ) perturbatively 



lmt 2 (s) = 0, lmt 4 (s) = a(s)\t 2 (s)\ 2 , ... => lm^- = a(s), (2.2) 



and cannot generate poles. Therefore the resonance region lies beyond the reach of standard ChPT. 
This region however, can be reached combining ChPT with dispersion theory either for the ampli- 
tude [|6|] or for the inverse amplitude through the IAM [J7|, ||, [J. 

Here we review our recent derivation [^] of a modified IAM formula to properly account for 
the Adler zero region. The main ingredient is a dispersion relation for the inverse amplitude, whose 
analytic structure consist on a right cut from s = Am 2 K to s = °°, a left cut from s = — oo to s = 0, and 
possible poles comming from zeros of the amplitude. Indeed, for the scalar waves the amplitude 
vanishes at the so called Adler zero, sa , that lies in the real axis below threshold, thus within the 

ChPT region of applicability. Its position can be approximated from ChPT, i.e., sa = s 2 + S4 H , 

where t% vanishes at s 2 , t% + 1\ at s 2 + S4, etc. 

We now write a once substracted dispersion relation for the inverse amplitude, where the 
substraction point has been chosen to be the Adler zero, sa- 

1 d/ n.i/,(0 ^ (1/t)+PC(1/()| (23) 

t{s) K Jrc (s-sa){s-s) 

where LC stands for a similar integral over the left cut and PC stands for the contribution of the 
pole at the Adler zero, which reads: 



PC (1/0 = -(s-s A ) Res 



W) \ 1 t"( SA ) 



s >= SA \(s'-s)(s'-s A )J t'(s A )(s-s A ) It 1 (s A ) 



2 • 



The different terms in Eq. ( |2.3[ ), can be evaluated in the following way: 

• The right cut can be exactly evaluated taking into account the elastic unitarity condition 
Eqs. 0, Q, lm\/t(s') = -o(s') = -lmt 4 (s')/t 2 (s'), tors' G (4m|,°o). 

• The pole contribution only involves amplitude derivatives evaluated at the Adler zero, which 
is a low energy point, so they can be can be safely approximated with ChPT. 

• The left cut, which is 1 / (s' — s) suppressed for s values near the physical region, is weighted 
at low energies, so it is appropiate to approximate it with ChPT. 

• Finally, in the cut contributions we approximate (s — s A )/(s' — Sa) — (s — S2)/(s ! — s 2 ), which 
is its LO chiral expansion and a remarkably good approximation as long as s' is is sufficiently 
far from s 2 and sa, which is indeed the case for the cut integrals. 



Altogether we are able to write Eq. ( |2.3[ ) as 

-J- ~ -RC(t A /t 2 2 ) -LC(t A /t 2 2 )+PC{l/t) = -4p- +PC{t A /t 2 )+PC(\/t), (2.4) 
t(s) v ' q(s) 

-U(s)/4(s)+PC(t 4 /t%) 
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where we have taken into account that a once subtracted dispersion relation can be written for U/t^, 
the substraction point being s 2 : 

4H = / ds 1 > \ > +LCitA tl) +PC(k &, (2.5) 

frjO*) 71 JRC [s' -s 2 )(s' -s) 



where PC{t^/t 2 ) stands for the contribution of the triple pole at s 2 . From Eq. (2^) one easily 
arrives at the modified IAM (mlAM) formula [§] : 

f™ 1 ®* -^ tl{ °\ r-rrr , A( S )= k (s 2 )- {S2 - SA){S - S2) (t 2 (s 2 )-t' 4 (s 2 )). (2.6) 



The standard IAM formula is recovered for A(s) = 0, which holds exactly for all partial waves 
except the scalar ones. In the original IAM derivation ||] A (s) was neglected, since it formally 
yields a NNLO contribution and is numerically very small, except near the Adler zero, where it 
diverges. However, if A(s) is neglected, the IAM Adler zero occurs at s 2 , correct only to LO, it is 
a double zero instead of a simple one, and a spurious pole of the amplitude appears close to the 
Adler zero on both the first and second sheet. All of these caveats are removed with the mlAM, Eq. 



(J2.6|). The differences in the physical and resonance region between the IAM and the mlAM are 
less than 1%. However, as we will see, for large m n the a poles appear as two virtual poles on the 
second sheet below threshold, one of them moving towards zero and approaching the Adler zero 
region, where the IAM fails. Thus, we will use for our calculations the mlAM, although it is only 
relevant for the mentioned second a pole, and only when it is very close to the spurious pole. 

In summary, there are no model dependences in the approach, but just approximations to a 
given order in ChPT Let us remark that, in our derivation, ChPT has been used always at low 
energies to evaluate parts of a dispersion relation, whose elastic unitarity cut is taken into account 
exactly. Thus, the IAM formula is reliable up to energies where inelasticities become important, 
even though ChPT does not converge at those energies, because ChPT is not being used there. This 
argument also holds when the a moves near threshold, because the Adler zero, which is the point 
where ChPT is used, is still at low energies far away from the pole. Note that the virtual pole that 
goes down in energies is on the second sheet, so does not prevent the use ChPT at the Adler zero, on 



the first sheet. Remarkably, the simple formula of the elastic mlAM, Eq.(2.6), (or the standard IAM 
one), while reproducing the ChPT expansion at low energies, is also able to generate both the p 
and a resonances with values of the LECs compatible with standard ChPT [^]. In other words, the 
IAM generates the poles ^ associated to these resonances in the second Riemann sheet. The fact 
that resonances are not introduced by hand but generated from first principles and data, is relevant 
because the existence and nature of scalar resonances is the subject of a long-lasting intense debate. 

To be precise, the IAM, when reexpanded, reproduces the ChPT series up to the order to which 
the input amplitude was evaluated and, in particular, the quark mass dependence agrees with that of 
ChPT up to that order. A few of the higher order terms are produced correctly by the unitarization 
but not the complete series — for a discussion of this issue for the scalar pion form factor see 
Ref.JT0|]. However, the formalism just described still provides us with a fair estimate of the quark 
mass dependence of the resonance properties. 
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Figure 1: Left: Movement of the a (dashed lines) and p (dotted lines) poles for increasing m n (direction 
indicated by the arrows) on the second sheet. The filled (open) boxes denote the pole positions for the a (p) 



at pion masses m n = 1,2, and 3 x m^ hys , respectively. For m n = 3m v ^ ly " three poles accumulate in the plot 
very near the n% threshold. Note that all poles are always far enough from the Adler zero (circles). Right: 
Comparison of our results for the M p dependence on m n with some recent lattice results from [flT|], The grey 
band covers only the error coming from the LECs uncertainities. 



3. Results 

By changing m n in the amplitudes we see how the poles generated with the IAM evolve. We 
will use the LECs values 10 3 / 3 r = 0.8 ± 3.8 and 10 3 /£ = 6.2 ± 5.7 from [JTJ] and fit the mlAM to 
data up to the resonance region to find 10 3 /[ = —3.7 ±0.2 and lO 3 ^ = 5.0 ±0.4. These LECs are 
evaluated at \x = 770 MeV. 

The values of m n considered should fall within the ChPT range of applicability and allow for 
some elastic %% regime below KK threshold. Both criteria are satisfied if m K < 500 MeV, since 
SU (3) ChPT still works with such kaon masses, and because for m K ~ 500 MeV, the kaon mass 
becomes ~ 600, leaving 200 MeV of elastic region. Of course, we expect higher order corrections, 
which are not considered here, to become more relevant as m % is increased. Thus, our results 
become less reliable as m n increases due to the 0{p^) corrections which we have neglected 

Fig. [I] (left) shows the evolution of the a and p pole positions as m n is increased. In order to 
see the pole movements relative to the two pion threshold, which is also increasing, all quantities 
are given in units of m n , so the threshold is fixed at ^fs = 2. Both poles moves closer to threshold 
and they approach the real axis. The p poles reach the real axis as the same time that they cross 
threshold. One of them jumps into the first sheet and stays below threshold in the real axis as a 
bound state, while its conjugate partner remains on the second sheet practically at the very same 
position as the one in the first. In contrast, the a poles go below threshold with a finite imaginary 
part before they meet in the real axis, still on the second sheet, becoming virtual states. As m K is 
increased further, one of the poles moves toward threshold and jumps through the branch point to 
the first sheet and stays in the real axis below threshold, very close to it as m n keeps growing. The 
other a pole moves down in energies further from threshold and remains on the second sheet. This 
analytic structure, with two very asymmetric poles in different sheets for a scalar wave, could be 
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Figure 2: m % dependence of resonance masses (left) and widths (right) in units of the physical values. 
In both panels the dark (light) band shows the results for the a (p). The width of the bands reflects the 
uncertainties induced only from the uncertainties in the LECs. The dotted line shows the a mass dependence 
estimated in Ref. Jlq]. The dashed (continuous) line shows the m K dependence of the a (p) width from the 
change of phase space only, assuming a constant coupling of the resonance to nn. 



a signal of a prominent molecular component [12, 13], at least for large pion masses. Similar pole 
movements have been also found within quark models [14] and in finite density analysis [15]. 

In Fig. @ (left) we show the m n dependence of M p and M a (defined from the pole position 
^/Spoie = M— iT/2), normalized to their physical values. The bands cover the LECs uncertainties. 
We see that both masses grow with m n , but M a grows faster than M p . Below m n ~ 330 MeV we 
only show one line because the two conjugate a poles have the same mass. Above 330 MeV, these 
two poles lie on the real axis with two different masses. The heavier pole goes towards threshold 
and around m % ~ 465 MeV moves into the first sheet. Note also that the m n dependence of M a is 
much softer than suggested in [16], shown as the dotted line, which in addition does not show the 
two virtual poles. 

In the right panel of Fig. |2| we show the m n dependence of T p and Y a normalized to their 
physical values, where we see that both widths become smaller. We compare this decrease with the 
expected reduction from phase space as the resonances approach the %% threshold. We find that 
Tp follows very well this expected behavior, which implies that the pun coupling is almost m K 
independent. In contrast r CT shows a different behavior from the phase space reduction expecta- 
tion. This suggest a strong m % dependence of the a coupling to two pions, necessarily present for 
molecular states [JT3], O]- 

Fig. [j] (right) shows our results for the p mass (here defined as the point where the phase shift 
crosses n/2, except for those m % values where the p becomes a bound state, where it is defined then 
from the pole position) dependence on m K compared with some recent lattice results [jll|], where 
we also quote the PDG value for the p mass. Taking into account the incompatibilities within errors 
between different lattice colaborations, we find a qualitaive good agreement with the lattice resuts. 
Also, we have to take into account that the m n dependence in our approach is correct only up to 
NLO in ChPT, and we expect higher order corrections to be important for large pion masses. 
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4. Summary 

We have reviewed our derivation of a modified version of the IAM [J20, which is derived from 
the first principles of analyticity, unitarity, and ChPT at low energies and has the correct behaviour 
around the Adler zero. It is able to generate the a and p resonance poles without any a priori 
assumptions, and yields the correct dependence on the pion mass up to NLO in ChPT. We have 
predicted the evolution of the resonance pole positions with increasing pion mass [J3j] and have seen 
that both resonances become bound states. We have also shown that the p Tin coupling constant is 
almost m n independent and we have found a qualitative agreement with some lattice results for the 
p mass evolution with m % . These findings might be relevant for studies of the meson spectrum and 



form factors — see Ref. J18Q — on the lattice. 
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